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ABSTRACT 

Context. Argonium has recently been detected as a ubiquitous molecule in our Galaxy. Model calculations indicate that its abundance 
peaks at molecular fractions in the range of lO^'* to 10“^ and that the observed column densities require high values of the cosmic 
ray ionization rate. Therefore, this molecular cation may serve as an excellent tracer of the very diffuse interstellar medium (ISM), as 
well as an indicator of the cosmic ray ionization rate. 

Aims. We attempted to detect ArH^ in extragalactic sources to evaluate its diagnostic power as a tracer of the almost purely atomic 
ISM in distant galaxies. 

Methods. We obtained ALMA observations of a foreground galaxy at z = 0.89 in the direction of the lensed blazar PKS 1830-211. 
Results. Two isotopologs of argonium, ^*ArH+ and ^*ArH+, were detected in absorption along two different lines of sight toward 
PKS 1830-211, known as the SW and NE images of the background blazar. The argonium absorption is clearly enhanced on the 
more diffuse line of sight (NE) compared to other molecular species. The isotopic ratio ^®Ar/^*Ar is 3.46 ± 0.16 toward the SW image, 
i.e., significantly lower than the solar value of 5.5. 

Conclusions. Our results demonstrate the suitability of argonium as a tracer of the almost purely atomic, diffuse ISM in high-redshift 
sources. The evolution of the isotopic ratio with redshift may help to constrain nucleosynthetic scenarios in the early Universe. 

Key words, quasars: absorption lines - galaxies: ISM - astrochemistry - galaxies: abundances - nucleosynthesis - quasars: individ¬ 
ual: PKS 1830-211 


1. Introduction 

Our knowledge of interstellar hydride species has been consid¬ 
erably increased by recent submillimeter missions, in particu¬ 
lar the Heterodyne Instrument for the Far-Infrared (HIFI) on 
board Herschel and the German REceiver At Terahertz frequen¬ 
cies (GREAT) on board the Stratospheric Observatory for Far- 
Infrared Astronomy (SOFIA), but also ground-based observa¬ 
tories, such as the Atacama Pathfinder Experi ment (APEX). 
The newly d etected hydride specie s are SH ^ (iMe nten et al 
201 ll) . OH^ (|W^rowsG_eL^ l20I0l) . HoO^ dOs senkonf et al. 


, H7CL (jnr et alTl20IOl) . HCL (iDe Luca et al.ll2012l). SH 


(iNeufeld et al.l 201^. and . most recently, ArH^ dBarlow et al.l 
l20I3t ISchilke et al.ll20I4ll . OH^ and H^O^ were also observed 
extensively in extragalactic sources (Ivan der Werf et ^ 120101: 

IWeiB et alJl2010t [Gonzalez-Alfonso et al. 11201^ 

The large column densities of OH^ with respect to H20^ 
were explained by both molecules residing preferentially in 
the largely atomic, diffuse interstellar medium (ISM), with the 
largest fractional abundances at a molecular fraction of around 
0.04, and their unexpectedly large column densities require cos¬ 
mic r ay ionization rates C considerably larger tha n in the dense 
ISM dHollenbach et al.ll2()12l:llndriolo et aLkoIl . 


The ArH^ cation was initially identified toward the Crab 
Nebula through its 7 = 1-0 and 2-1 transitions, with 
OH^ 1 - 0 as the only other emission feature in the spec¬ 
trum recorded with SPIRE ( Spectral and Photom etric Image RE- 
ceiver) on board He rschel dBarlow et al.ll2013l) . More recently, 
ISchilke et al.l d2014l) reported on their early detections of the 
ArH^ 7=1-0 transition in absorption toward six bright con¬ 
tinuum sources using Herschel-Y\\¥l. The absorption patterns 
were unique in their appearances; ArH^ showed up in all ve¬ 
locity components associated with diffuse foreground molecu¬ 
lar clouds, but was conspicuously absent at velocities related to 
the sources themselves. The observations were reproduced by 
models in which argonium is present only in the almost purely 
atomic diffuse ISM, with peak abundance at molecular frac¬ 
tions in the 10“^ to 10“^ range. Basically, high values of ^ fa¬ 
vor its formation, consistent with analyses of OH^ and H20^, 
which trace slightly higher, but still low H2 fractions (see, e.g., 
iHollenbach et al.ir2012illndriolo et aO20L5 ). whereas molecular 


fractions > 10“^ favor its destruction by the reaction with H2 to 
produce Ar and . Although ArH^ abhors molecular clouds, 
it does rely upon a small amount of H2 for its existence. More¬ 
over, with its specificity to the ISM with a molecular fraction of 
< 10'^, it is a much better tracer of the almost purely atomic ISM 
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than the 21 cm hyperfine structure line of atomic hydrogen, be¬ 
cause the latter is also seen in th e ISM with a molec ular fraction 
of 0.01 or even larger than 0.5 (ISchilke et al.ll201^ . Argonium 
may also be used to infer ^ in combination with other tracers. 


Because argonium is ubiquitous in the Galactic diffuse ISM, 
the next step is to search for it in extragalactic sources. The J - 
1 - 0 transition of36ArH+ at 617525.23 + 0.15 MHtQ is close to 
a strong atmospheric water line (532 -4-4i near 620.7 GHz0, and 
the ^*ArH^ transition is less than 900 MHz lower in frequency. 
Therefore, ground-based observations of ArH^ in absorption in 
nearby galaxies are very difficult. As Herschel is no longer in 
operation and SOFIA has no receiver at these frequencies yet, 
searches toward galaxies at low to moderate redshifts are more 
promising, given the current instrumental capabilities. 

The unnamed foregro und galaxy at z = 0.88582 

dWiklind & Combes! Il996h in the direction of the blazar 
PKS 1830-211 is particularly well suited for such searches. 
The galaxy is lensi ng the background radiation com¬ 
ing from the blazar (jfaamesh Rao & Subrahmany^ 1 19881: 
ISubrahmanvan et al.lll990ll : two images, to the SW and NE of 
the blazar and separated by 1", are especially strong. Absorp¬ 
tion spectra toward the SW and NE images probe two separate 
lines of sight across this foreground galaxy located on opposite 
sides of the nucleus. Eor most species, the absorption toward 
the SW image is much stronger than that toward the NE im¬ 
age, since it probes denser gas, with higher H 2 column densi¬ 
ties. The blazar itself appears to be devoid of line-absorption or 
emission at radio frequencies. The line of sight toward the SW 
image of PKS 1830-211 is, to date, the extragalacti c object with 
the largest nu mber of molecules detected (~40, see lMuller et al.l 
l2nilLl2ni4hh . 


The redshift z = 0.89 corresponds to a look-back time of 
~7.5 Gyr, when the Universe was slightly less than half its cur¬ 
rent age. At that time, the enrichment of the ISM with heavy 
elements was dominated by massive stars. This could lead to 
isotopic ratios of the elements different from those in the local 
galactic ISM. Molecular species may be much more promising 
in this regard than atoms, because of their more distinct spectro¬ 
scopic patterns, e.g. different rotational patterns, caused by dif¬ 
ferences in the reduced masses. Indeed, isotopic ratios have been 
determined in the foreground galaxy for molecules containing 
moderately heavy elements such as Si, S, and Cl in addition to 
C, N, a nd O, and show some differences compared to their solar 
values (iMuller et al.ll2005l201 lll2014al) . 


The Atacama Large Millimeter/sub millimeter Array 
(ALMA) now offers high sensitivity and wide frequency cov¬ 
erage for spectroscopic studies. Already in the Early Science 
Cycle 0, despite the limited number of antennas and only four 
spectral tunings, a large number (17) of common interstellar 
species, such as CO, CH, H 2 O, HCO+, HCN, C 2 H, and NH 3 , 
were observed toward PKS 1 830-211, including two new 
extra galactic dete ctions, H 2 CU (iMuller et al.ll20143) and NH 2 
(see Muller et~^ l2014bl for the presentation of the ALMA 
Cycle 0 survey and overall results). 


The present Letter reports the first extragalactic detec¬ 
tion of ^®ArH^ and ^®ArH^ from ALMA observations toward 
PKS 1830-211. 




Velocity (km/s, heliocentric frame, z=0.88582) 

Fig. 1. Spectra of 36ArH+ and 3*ArH+ toward PKS 1830-211 SW (top) 
and NE (bottom) images, together with two para-fliO* fine structure 
lines of the To - loi transition; none of the transitions shows hyperfine 
splitting. 


2. Observations and spectroscopic data 

Observations reported here were obtained with ALMA on 2015 
May 19 (Cycle 2), in about 20 min on-source time under 
very good atmospheric conditions (precipitable amount of wa¬ 
ter vapor of ~0.3 mm) and with receivers tuned to ~327 GHz 
(Band 7fl The three ArH^ isotopologs were observed simul¬ 
taneously in the same 1.875 GHz wide spectral window. Two 
additional spectral windows were centered at ~337 GHz to 
cover the HzO^ lio - loi, J - 0.5 - 1.5 line (rest frequency 
634.27 GHz) and at ~339 GHz (no lines were detected in this 
window, which was used to map the continuum). The bandpass 
response of the antennas was calibrated using the radio-bright 
quasar J1924-2914, and the gain solutions versus time were de¬ 
termined from observations of the quasar J1832-2039, within 
1° on the sky from PKS 1830-211. The overall standard cali¬ 
bration was applied by the ALMA pipeline. One more step of 
self-calibration was performed and calibrated visibilities w ere 
then fit with the task UVMULTILIT (iMartf- Vidal et akllMl to 
extract the spectra toward the point-like SW and NE images of 


’ From the CDMS catalog (IMuller et all 1200 IL l2005h . see also sec- 
tion|2] 

^ See, e.g., the JPL catalog dPickett et al.lll998h . 

^ At z = 0.89, the frequencies of the 7 = 1-0 line of ArH^ iso¬ 
topologs are redshifted close to the Sij - 422 atmospheric water line 
near 325.2 GHi^l. 
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H. S. P. Muller et al.: Detection of extragalactic argonium, ArH^, toward PKS 1830-211 
Table 1. Spectroscopic data" of the y = 1 - 0 transitions of ArH+ isotopologs observed in the present study. 


Species 

Rest frequency 
(MHz) 

A ■ 

(10-3 s-') 

^low 

(K) 

JrdviSWy 

(kms-') 

Av(SW)‘^ 

(kms-') 

Voff(SW)^ 

(kms-') 

JrdviSWy 

(kms-') 

jT£/v(NE)« 

(kms"') 

'''’ArH+ 

617525.23 (15) 

4.33 

0.0 

13.65 (17) 

57.3 (8) 

5.6 (4) 

15.88 (21) 

7.61 (12) 

38ArH+ 

616648.76 (8) 

4.31 

0.0 

3.90(15) 

57.3(8) 

5.6 (4) 

4.59 (21) 

1.68(12) 

40ArH+ 

615858.15(5) 

4.30 

0.0 




< 0.42'’ 



Notes. “ Numbers in parentheses are one standard deviation in units of the least significant figures. * Einstein A value. ^ Velocity integrated optical 
depth toward the SW image from a single Gaussian fit. Full width at half maximum of main SW image feature from that fit constrained to be 
equal for all three species. Offset of main SW image centroid velocity with respect to the systemic velocity from that fit constrained to be equal 
for all three species, f Optical depth toward the SW image integrated over channels from -80 to +110 kms“*. ® Optical depth toward the NE 
image integrated over channels from -280 to -80 km s“*. * Upper limit, 3cr uncertainty. 


PKS 1830-211. The two images (separated by ~1") were well 
resolved by the synthesized beam of ~0.5". The velocity reso¬ 
lution is 1.0 kms^' after Hanning smoothing. Two atmospheric 
lines due to O3 (at 326.901 and 327.845 GHzj^show up in the 
spectra (Fig.[T]l. They are easily identified and disentangled from 
absorption intrinsic to the source, since they occur at the same 
frequency in both S W and NE spectra, see Fig. lA.ll 

Spec troscopic data of ArH^ tr ansitions were taken from the 
CDMS dMiiller et al.ll200lL l2005h and are summarized in Ta¬ 
ble [T] An isotopic invariant fit of div erse rotational and rovib ra- 
tional data was carried out, similar to lOdashima et al.l (Il999ll . to 
determine spectroscopic parameters of the molecule and to de¬ 
rive rest frequencies of rotational transitions of argonium iso¬ 
topologs. The most relevant data for our observations are the 
J = 1-0 transition s of ^®ArD^, ^*ArD^, and "^‘’ArD^ from 
iBowman et al.l (Il983h . As "^°Ar is the most abundan t Ar iso¬ 
tope on Earth, almost all of the remaining rotatio nal (iLiu et al.l 
Il987t iBrownetal. 1988t Odashima et al J l999h and rovibra- 


tional data ( Brault & Davis 19821 : |johnsll984 Filgueira & BlomI 


Il988ll used in that calculation refer to "^°ArH+ or to "^°ArD 


The "^°ArH^ ground state dipole moment is 2.18 D (ICheng et al.l 
I2007h . 


3. Results and discussion 

3.1. Column densities and abundances 

Absorption signatures due to ^®ArH^ and ^^ArH^ were detected 
with good to very good signal-to-noise ratios toward both im¬ 
ages of PKS 1830-211, SW and NE, see FiglT] The images and 
overview spectra are shown in Fig. lA.ll in the Appendix. 

Toward the SW image, a simple fit with one Gaussian veloc¬ 
ity component and the same width and velocity offset for both 
isotopologs yields a centroid velocity of 5.6 + 0.4 kms * and a 
full width at half maximum of 57.3 + 0.8 km s ', which, to first 
order, reproduces well both ^®ArH^ and ^^ArH^ absorptions. In¬ 
tegrated opacities of 13.64 + 0.19 and 3.89 + 0.16 kms * were 
derived for ^®ArH^ and ^^ArH^, respectively. We obtain only an 
upper limit of 0.42 kms ' (at 3cr confidence level) in the case 
of ""'ArH^. However, the ArH^ line profile shows a slight devi¬ 
ation from a simple Gaussian, and fit results and values of inte¬ 
grated opacity over velocity channels from -80 to +110 km s“* 
are given in a separate column in Table [T] The weak velocity 
component at +170 kms ', previously only detected in the lines 
of HCO'^, H CN, and H^O and all three strongly saturated near 
V = 0 km s“' (iMuller et al.l201 lLl2014bl) . is - surprisingly - also 
detected in the ArH^ spectrum. 


In contrast to the SW line of sight, where the absorption 
mostly resides in one bulky absorption feature (with the excep¬ 
tion of the additional +170 kms * component), the ArH^ ab¬ 
sorption profile toward the NE image shows a remarkable series 
of narrow (a few kms^') features spanning over ~200 kms~', 
also se en, e.g., in the absorption profile of H2O (iMuller et al.l 
I2014hh . 

With the assumptions that the rotational excitation temper¬ 
ature is equal to the radiation temperature of the cosmic mi¬ 
crowave background at z = 0.89, 5.14 K, and that the lines are 
optically thin, the total column densities of ^®ArH^ along the SW 
and NE lines of sight are 2.7 x lO'^ cm“^ and 1.3 x 10'^ cm~^ 
respect ively, i.e., differing by only a factor ~2. IMuller et al.l 
(l2014ah found a similar ratio of ~3 for H2CE. In contrast, the 
H2 column density, estimated from proxies such as CH, is abou t 
one order of magnitude higher toward SW (IMuller et al.ll2014'5) . 
The NE lin e of sight is known to be mo re diffuse and richer in 
atomic gas dKoopmans & de Bruvnl200^ . thus the enhancement 
of ArH^ on the NE line of sight relative to the S W one is not sur¬ 
prising. The detection of the ArH^ absorption in the +170 km s“* 
velocity component toward the SW image suggests that this re- 
gion has a low molecula r fraction as well, or that time variations 
dMuller & Guelinll2008l) have made this feature stronger. 

Since it is expected that ArH^ traces a gas component with 
a very low H2/H fraction, it would seem natural to compare 
its absorption profile with that of Hi. The H i absorption has 
been previously observed toward t he blazaiQ dChengalur et al.l 
ll999l : f^oDmans & de Bruvnl2005h . but only at low angular res¬ 
olution, not resolving the background continuum morphology. 
PKS 1830-211 is radically different at centimeter and submm 
wavelengths: the NE and SW images, which are point-like in 
the submm, are more extended in the cm: in addition, the emis- 
sion from the pseudo Einste in ring can be seen in the cm regime 
dSubrahmanvan et^Il990l) . Thus, comparison of spectra be¬ 
tween such different frequencies is not straightforward. 

Both H20 ' ^ and ArH^ trace the diffuse ISM component 
dindriolo et alJ l2015l) . Therefore, one may assume that both 
species show similar absorption features. This is not the case 
toward either of the PKS 1830-211 images, as can be seen in 
Fig.m this also demonstrated in a correlation plot (Fig. IA.2b in 
the Appendix. The 634 GHz line of H20'^ was observed simul¬ 
taneously with ArH^, while the 607 GHz line, with a higher rel¬ 
ative intensity and better signal-to-noise ratio, was observed in 
2014 (Muller et al., in preparation). There appear to be no large 
variations in the H20^ spectra, which were observed about one 


See also http://www.atnf.csiro.au/projects/askap/news_commissioning 
_10042014.jpg 
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year apart. A poor correlation between the two mo lecular ions 
is also seen in Galactic sources (ISchilke et al.ll2014l) . suggesting 
that the diffuse ISM with a molecular fraction 10“^ to 10“^ (as 
traced by ArH^) is different from the diffuse ISM with a molec¬ 
ular fraction around 0.04 (as traced by H20^). It is conceivable 
that at least part of ArH^ absorption arises in the warm (8000 K) 
ionized medium (WIM), rather than entirely in the least molecu¬ 
lar parts of the diffuse, neutral ISM. A detailed discussion of the 
chemistry of the WIM will be presented elsewhere. 


3.2. The ^^Ar/^^Ar isotopic ratio atz - 0.89 

The ^®Ar/^*Ar ratio is 3.46 + 0.16 (integrated over velocities 
from -80 to -i-llO kms“'; 3.50 + 0.14 from the Gaussian 
fit) toward the S W image, significa ntly lower than the solar 
value 5.50 + 0. 01 (IVogel et alj]201 ih, or the terrestrial value of 
5.305 + 0.068 (iBerglund & Wiesedlioi ih . Although ^°Ar (from 
the radio-active decay of is dominant on Earth, we note 
that it plays only a minor role in the ISM. The a elements sul¬ 
fur and silicon also sh ow ratios about two times lower than their 
terrestrial values (see iMuller et al.ll2006l |20 1 lb . In contrast, the 
^^Cl/^^Cl ratio was f ound to be identical to the terrestrial value 
(iMuller et al.ll2014il) . The ^®Ar/^*Ar ratio toward the NE image 
is 4.53 + 0.33, between the SW ratio and the solar and terrestrial 
values and compatible with either within the uncertainties. The 
3cr limit for the ^^ArH^/"^‘^ArH^ rati o of > 11 is not co nstraining 
considering the solar ratio of ~615 (I Vogel et al.ll201 Ih . which is 
ass umed to be close to the value in the local ISM. 

iKobavashi et al.l (l2006l) calculated nucleosynthesis yields of 
core-collapse supernovae (SNe) and hypernovae (HNe, which 
are more energetic by an order of magnitude) to predict the evo¬ 
lution of isotopic ratios in the Milky Way, SNe and HNe being 
the dominant contributors to elements from Na to Ee with the 
possible exception of a few selected isotopes. They presented 
yields for different progenitor masses and metallicities, and sub¬ 
sequently incorpor ated the results into mo dels for the enrichment 
of the Milky Way (IKobavashi et al.ll201 ll) . While the nucleosyn- 
thetic yields for many elements, including Ar and Ar, show 
a non-monotonic dependence on progenitor mass, the predicted 
^®Ar/^^Ar ratio tends to decline with increasing metallicity and 
is typically smaller for HNe than for SNe. Nevertheless, even 
for the highest metallicities considered in the Milky Way en¬ 
richment models (Z = 0.02, or solar metallicity), the predicted 
Ar/^®Ar ratio is larger than the observed ratio in either the Sun 
or the PKS 1830-211 absorber. There are several uncertain as¬ 
sumptions in the nucleosynthesis and Galactic enrichment mod¬ 
els, including the importance of convective mixing in the pro¬ 
genitor stars and the relative frequency of HNe; observations of 
isotopic ratios in diverse environments, such as those presented 
here, promise to provide valuable constraints for future mod¬ 
els. There is no evidence for low metallicity in the foreground 
galaxy, but the exact value cannot be determined at present. 


4. Conclusion 

We have detected for the first time extragalactic argonium in both 
of its important interstellar isotopologs. Therefore, ArH^ is not 
only a good tracer of the almost purely atomic diffuse ISM in 
the Milky Way, but it is also suitable for investigations of extra- 
galactic sources, possibly to a similar extent as OH^ or H20^. It 
will be interesting to search for ArH^ in a selection of galaxies 
at moderately low to higher redshifts to obtain a more complete 
picture of the atomic component of the extragalactic ISM and to 
look for possible trends in the Ar/^*Ar ratio with redshift. The 


Si, S, Cl, and Ar isotopic ratios appear to be promising probes 
into the nucleosynthesis in the early Universe because their yield 
ratios depend strongly on the type of supernova, and all four el¬ 
ements form molecules that have been detected even in distant 
galaxies. It will be interesting to see if these ratios correlate with 
f as both aspects are affected by SN activity. Searches for argo¬ 
nium in nearby galaxies may be more promising with the poten¬ 
tial future GREAT receiver on board SOEIA. Combined obser¬ 
vations of ArH^ with OH^, H20^, and HE or CH will provide in¬ 
sight into the partitioning of atomic and molecular hydrogen and 
into the cosmic ray ionization rate. Such an analysis is currently 
under way for PKS 1830-211 and will be published separately. 
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Appendix A: Complementary figures 
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Fig. A.l. Map of the 339 GHz continuum emission of PKS 1830-211 showing the two resolved lensed images of the blazar (left). Overview of 
the ArH+ absorption spectrum toward PKS 1830-211 SW (bottom right) and NE (top right) from the current observations. Expected positions of 
the main absorption feature are indicated for ^®ArH^, ^*ArH^, and ''°ArH^; the last is not detected. The noise increases toward lower frequencies 
because of the proximity of the atmospheric water line near 325.2 GHz. 



N(ArH+) [normalized] 


Fig. A.2. Normalized correlation plot of H20^ versus ^®ArH+ toward PKS 1830-211 SW, NE, and Sagittarius B2(M). The column densities 
of the two species have been normalized by the maximum values for each component and for each velocity channel. A perfect correlation (i.e. 
just a constant scaling factor between the two species) would result in clustering along the yellow line. The distributions of the two cations are 
poorly correlated toward the sources. The PKS 1830-211 data are from this work a nd from Muller et al., in preparation for ^®ArH^ and H20^, 
respectively; the corresponding Sagittarius B2(M) data are from ISchilke et al.l l l2014h and lSchilke et al.l ( l2013h . 
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